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Abstract 

The eigenvalue bounds obtained earlier [J. Phys. A: Math. Gen. 31 (1998) 963] 
for smooth transformations of the form V(x) = g(x 2 ) + f(—s) are extended to N- 
dimensions. In particular a simple formula is derived which bounds the eigenvalues 
for the spiked harmonic oscillator potential V(x) = x 2 + a > 0, A > 0, and is valid 
for all discrete eigenvalues, arbitrary angular momentum I and spatial dimension N. 
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Recently, a simple formula that bounds the eigenvalues E n of Schor- dinger's equation 

-^" + (g( x 2)+f(\))i; = E n i;, V(0) = 0, (1) 

where g and / are two smooth transformations of x 2 and respectively, was obtained 
by the present authors. They showed that E n can be approximated by the expression 



E n ^{ g ^)-^) + f^)-^) + 
Vsf) (4n + 2 + ^ 4 /'(l) + l) | , 



(2) 

n = 0,1,2,... 



This formula provides a lower bound = >) or an upper bound = <) to the 
exact eigenvalues i? n of Eq.(l) according as the transformation functions g and / are 
both convex (g" > 0, /" > 0) or both concave (g" < 0, /" < 0). This allowed us, 
for example, to obtain simple expressions which bound the spectrum of the spiked 
harmonic oscillator potential V(x) = Xx 2 + a > 1, n = 0, 1, 2, . . . , namely 

E n « e n (t) = (1 - + 2Xt 2 + 2\f\~(2n + 1) (3) 

where t is the real positive root of 

2/jat 2 ~ a - 4At 4 + 1 = 0. 

Here e n (t) is lower bound to E n when a > 2 and an upper bound when a < 2. The 
purpose of the present paper is to extend these results to the iV-dimensional case with 
arbitrary angular momentum number /. 
We notice first that the exact eigenvalues [3] 

E nl = VX{An + 2 + y/Afji + (21 + l) 2 ), n = 0,1,2,... 

of Schrodinger's equation with the Gol'dman and Krivchenkov potential V(x) = Xx 2 + 
-^2 in 3-dimensions can be extended [4] to the N dimensional case by replacing / with 
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I + — §• Indeed, these exact solutions could be generated from the well known 
solutions of harmonic oscillator potential by two simple transformations: first replace 
the angular momentum / in the harmonic oscillator energy expression V\(4n + 21 + 
3), n = 0, 1, 2, . . . by -\ + yV + (I + ^) 2 ! then replace / with I + ^ - §. Thus the 
exact eigenvalues of the iV-dimensional schrodinger equation with the Gol'dman and 
Krivchenkov potential are 

E { n f } = 2^(2?! + 1 + v/^ + (/ + A^/2-l) 2 ), n, I = 0, 1, 2, . . . (4) 

The method used to develop the results Eq.(5) and Eq.(6) of Ref.[l] can now be 
followed, but instead of formula (2) of Ref.[l] we use Eq.(4) above and obtain 

« miners,*) 

s,i>U 

where 



+ 2v^)(2n + 1 + y//'(l) + (/ + iV/2-l)2) J . 
The case where g(x 2 ) = \x@ and f{-^) = implies, from (5), that 



(5) 



^(M) = A(l-§)-* + (l-f)£ 



+ VSA/^" 2 + 1 + yj + (/ + iV/2 - 1 



(6) 



In particular, for the spiked harmonic oscillator potential with (3 = 2 it follows from 
(5) that the eigenvalue approximation is given by 

e^\t) = + 2At 2 + 2v/A(2n + 1) (7) 

where t is real positive root of 

2At 4 - ixat 2 ~ a - 2(1 + N/2 - l) 2 = 0. (8) 
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We now prove that for optimal t there is only one positive real root given by (8). If 
we let h(t) = 2At 4 - / uat 2 - Q -2(Z + iV/2-l) 2 , then for a < 2: h(t) -> -2(Z + iV/2-l) 2 
as t — > and — > oo as t — > oo. On the interval (0, oo) the function ft(£) has only 
one minimum occurring at 



Consequently for a < 2, Eq.(8) has only one real positive root. For a > 2, — > — oo 
as t — > + and to +oo as t — > oo. On the interval (0, oo), /i(£) is monotone increasing 
on the open interval (0, oo) and we conclude that (8) has only one real positive solu- 
tion for all a. The above discussion leads to the following simple expression for the 
energy bound approximations for the spiked harmonic oscillator potential valid for all 
dimensions N > 2, arbitrary angular momentum / > 0, and n > 



In Table 1 we exhibit the upper bounds Eq obtained by use of formula (7) for di- 
mensions N = 2 to 10 with a = 1.9, A = 1, and fx = 10, along with some accurate 
values obtained by direct numerical integration of Schrodinger's equation. Similar 
accurate numerical results could also be obtained by the use of perturbation methods 
such as the renormalized hypervirial perturbation method of Killingbeck [5]. In Table 
2 we exhibit the corresponding lower bounds obtained by use of formula (7) for 
dimensions iV = 2 to 10 with a = 2.1, A = 1, and n = 10. 

For the particular test problem discussed here, other approximation methods might 
also be considered. For example, if we let E(a) represent an eigenvalue of the operator 
H(a) = —A + x 2 + [ix~ a with \x fixed, then an immediate first-order approximation 
'formula' is provided by 





(9) 



E(a) w E(2) + (a-2)E'{2). 



(10) 



Eigenvalue bounds 



Page 5 



The problem now is to find E'(2). Since E(a) = (ip(a), H(a)t(j(a)), differentiation 
with respect to a and the minimal property of the expectation (i(^(a), H(2)ifj(a)) with 
respect to a leads to the expression: 

E'(2) = -M^(2),log0r)*- 2 V>(2)). (11) 

As an illustration of this result we consider the first and last lines of Table 1. We 
find for iV = 2 that E'{2) « -1.557 and £(1.9) = 8.4803; meanwhile for N = 10, 
E'(2) « -1.498 and £(1.9) w 12.3479. The same reasoning and method can be 
applied to Table 2. However, these results are particular to the example chosen for 
illustration, and they do not in general provide energy bounds. Given the correct 
convexity of the transformation functions, the geometrical methods described in the 
present paper provide energy bounds on all the discrete eigenvalues in all dimensions 
N>2. 

By extending the scope to N dimensions we have generalized our simple general 
eigenvalue approximation formulae for the potential 

V(x) = g(x 2 ) + f(±), 
x z 

where g and / are smooth monotone transformations of x 2 and respectively. These 
results may be used for exploratory purposes and also for seeding direct numerical 
methods. In Figs.(l) and (2) we show the potential, the eigenvalues, and the unnor- 
malized radial wave functions corresponding to the data in Tables (1) and (2). The 
computation of such results is greatly helped by a priori knowledge of the approximate 
location of the eigenvalues. 
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Table(l): Upper bounds Eq using (7) for H = —A + x 2 + -£r& for dimension N = 
2 to 10. The 'exact' values E 00 were obtained by direct numerical integration of 
Schrodinger's equation. 



N E o E Q0 

2 8.485 38 8.511 90 

3 8.564 36 8.590 21 

4 8.795 44 8.819 47 

5 9.163 09 9.184 61 

6 9.646 70 9.665 48 

7 10.225 05 10.241 20 

8 10.879 08 10.892 89 

9 11.592 98 11.604 78 

10 12.354 18 12.364 29 
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Table(2): lower bounds E^\ using (7) for H = —A + x 2 + -^rr for dimension N = 
2 to 10. The 'exact' values E 2 \ were obtained by direct numerical integration of 
Schrodinger's equation. 



N 



E. 



21 



E21 



2 
3 
4 
5 
6 
7 
8 
9 

10 



16.457 73 
16.826 41 
17.312 54 
17.895 07 
18.554 81 
19.275 58 
20.044 44 
20.851 25 
21.688 22 



16.543 63 
16.904 44 
17.381 71 
17.955 44 
18.607 07 
19.320 69 
20.083 41 
20.885 02 
21.717 61 
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Fig. (2) Graph of the eigenvalues E = E^\ for the Schrddinger equation with the po- 
tential V(x) = x 2 + ^tt and corresponding unnormalized wavefunctions in dimensions 
N = 2 (bottom) to 10 (top). 



